%% distribution of omega
load('./data/FigS7.mat')
y=unique(omega_save(:)/1e6);
histogram(y)

h=gca();
xlabel('\omega_r (MHz)')
ylabel('occurences')

h.FontSize=18;

% saveas(gcf,'FigS7','epsc')
